Folding of prestin's anion-binding site and the mechanism of outer hair cell electromotility.

Prestin responds to transmembrane voltage fluctuations by changing its cross-sectional area, a process underlying the electromotility of outer hair cells and cochlear amplification. Prestin belongs to the SLC26 family of anion transporters yet is the only member capable of displaying electromotility. Prestin's voltage-dependent conformational changes are driven by the putative displacement of residue R399 and a set of sparse charged residues within the transmembrane domain, following the binding of a Cl- anion at a conserved binding site formed by the amino termini of the TM3 and TM10 helices. However, a major conundrum arises as to how an anion that binds in proximity to a positive charge (R399), can promote the voltage sensitivity of prestin. Using hydrogen-deuterium exchange mass spectrometry, we find that prestin displays an unstable anion-binding site, where folding of the amino termini of TM3 and TM10 is coupled to Cl- binding. This event shortens the TM3-TM10 electrostatic gap, thereby connecting the two helices, resulting in reduced cross-sectional area. These folding events upon anion binding are absent in SLC26A9, a non-electromotile transporter closely related to prestin. Dynamics of prestin embedded in a lipid bilayer closely match that in detergent micelle, except for a destabilized lipid-facing helix TM6 that is critical to prestin's mechanical expansion. We observe helix fraying at prestin's anion-binding site but cooperative unfolding of multiple lipid-facing helices, features that may promote prestin's fast electromechanical rearrangements. These results highlight a novel role of the folding equilibrium of the anion-binding site, and help define prestin's unique voltage-sensing mechanism and electromotility.

Folding of Prestin's Anion-Binding Site and the Mechanism of Outer Hair Cell Electromotility.

Prestin responds to transmembrane voltage fluctuations by changing its cross-sectional area, a process underlying the electromotility of outer hair cells and cochlear amplification. Prestin belongs to the SLC26 family of anion transporters yet is the only member capable of displaying electromotility. Prestin's voltage-dependent conformational changes are driven by the putative displacement of residue R399 and a set of sparse charged residues within the transmembrane domain, following the binding of a Cl - anion at a conserved binding site formed by amino termini of the TM3 and TM10 helices. However, a major conundrum arises as to how an anion that binds in proximity to a positive charge (R399), can promote the voltage sensitivity of prestin. Using hydrogen-deuterium exchange mass spectrometry, we find that prestin displays an unstable anion-binding site, where folding of the amino termini of TM3 and TM10 is coupled to Cl - binding. This event shortens the TM3-TM10 electrostatic gap, thereby connecting the two helices, resulting in reduced cross-sectional area. These folding events upon anion-binding are absent in SLC26A9, a non-electromotile transporter closely related to prestin. Dynamics of prestin embedded in a lipid bilayer closely match that in detergent micelle, except for a destabilized lipid-facing helix TM6 that is critical to prestin's mechanical expansion. We observe helix fraying at prestin's anion-binding site but cooperative unfolding of multiple lipid-facing helices, features that may promote prestin's fast electromechanical rearrangements. These results highlight a novel role of the folding equilibrium of the anion-binding site, and helps define prestin's unique voltage-sensing mechanism and electromotility.

Factors That Control the Force Needed to Unfold a Membrane Protein in Silico Depend on the Mode of Denaturation.

Single-molecule force spectroscopy methods, such as AFM and magnetic tweezers, have proved extremely beneficial in elucidating folding pathways for soluble and membrane proteins. To identify factors that determine the force rupture levels in force-induced membrane protein unfolding, we applied our near-atomic-level Upside molecular dynamics package to study the vertical and lateral pulling of bacteriorhodopsin (bR) and GlpG, respectively. With our algorithm, we were able to selectively alter the magnitudes of individual interaction terms and identify that, for vertical pulling, hydrogen bond strength had the strongest effect, whereas other non-bonded protein and membrane-protein interactions had only moderate influences, except for the extraction of the last helix where the membrane-protein interactions had a stronger influence. The up-down topology of the transmembrane helices caused helices to be pulled out as pairs. The rate-limiting rupture event often was the loss of H-bonds and the ejection of the first helix, which then propagated tension to the second helix, which rapidly exited the bilayer. The pulling of the charged linkers across the membrane had minimal influence, as did changing the bilayer thickness. For the lateral pulling of GlpG, the rate-limiting rupture corresponded to the separation of the helices within the membrane, with the H-bonds generally being broken only afterward. Beyond providing a detailed picture of the rupture events, our study emphasizes that the pulling mode greatly affects the factors that determine the forces needed to unfold a membrane protein.

Challenges and Advantages of Accounting for Backbone Flexibility in Prediction of Protein-Protein Complexes.

Predicting protein binding is a core problem of computational biophysics. That this objective can be partly achieved with some amount of success using docking algorithms based on rigid protein models is remarkable, although going further requires allowing for protein flexibility. However, accurately capturing the conformational changes upon binding remains an enduring challenge for docking algorithms. Here, we adapt our Upside folding model, where side chains are represented as multi-position beads, to explore how flexibility may impact predictions of protein-protein complexes. Specifically, the Upside model is used to investigate where backbone flexibility helps, which types of interactions are important, and what is the impact of coarse graining. These efforts also shed light on the relative challenges posed by folding and docking. After training the Upside energy function for docking, the model is competitive with the established all-atom methods. However, allowing for backbone flexibility during docking is generally detrimental, as the presence of comparatively minor (3-5 Å) deviations relative to the docked structure has a large negative effect on performance. While this issue appears to be inherent to current forcefield-guided flexible docking methods, systems involving the co-folding of flexible loops such as antibody-antigen complexes represent an interesting exception. In this case, binding is improved when backbone flexibility is allowed using the Upside model.

Prediction and Validation of a Protein's Free Energy Surface Using Hydrogen Exchange and (Importantly) Its Denaturant Dependence.

The denaturant dependence of hydrogen-deuterium exchange (HDX) is a powerful measurement to identify the breaking of individual H-bonds and map the free energy surface (FES) of a protein including the very rare states. Molecular dynamics (MD) can identify each partial unfolding event with atomic-level resolution. Hence, their combination provides a great opportunity to test the accuracy of simulations and to verify the interpretation of HDX data. For this comparison, we use Upside, our new and extremely fast MD package that is capable of folding proteins with an accuracy comparable to that of all-atom methods. The FESs of two naturally occurring and two designed proteins are so generated and compared to our NMR/HDX data. We find that Upside's accuracy is considerably improved upon modifying the energy function using a new machine-learning procedure that trains for proper protein behavior including realistic denatured states in addition to stable native states. The resulting increase in cooperativity is critical for replicating the HDX data and protein stability, indicating that we have properly encoded the underlying physiochemical interactions into an MD package. We did observe some mismatch, however, underscoring the ongoing challenges faced by simulations in calculating accurate FESs. Nevertheless, our ensembles can identify the properties of the fluctuations that lead to HDX, whether they be small-, medium-, or large-scale openings, and can speak to the breadth of the native ensemble that has been a matter of debate.

Lipid bilayer induces contraction of the denatured state ensemble of a helical-bundle membrane protein.

Defining the denatured state ensemble (DSE) and disordered proteins is essential to understanding folding, chaperone action, degradation, and translocation. As compared with water-soluble proteins, the DSE of membrane proteins is much less characterized. Here, we measure the DSE of the helical membrane protein GlpG of Escherichia coli (E. coli) in native-like lipid bilayers. The DSE was obtained using our steric trapping method, which couples denaturation of doubly biotinylated GlpG to binding of two streptavidin molecules. The helices and loops are probed using limited proteolysis and mass spectrometry, while the dimensions are determined using our paramagnetic biotin derivative and double electron-electron resonance spectroscopy. These data, along with our Upside simulations, identify the DSE as being highly dynamic, involving the topology changes and unfolding of some of the transmembrane (TM) helices. The DSE is expanded relative to the native state but only to 15 to 75% of the fully expanded condition. The degree of expansion depends on the local protein packing and the lipid composition. E. coli's lipid bilayer promotes the association of TM helices in the DSE and, probably in general, facilitates interhelical interactions. This tendency may be the outcome of a general lipophobic effect of proteins within the cell membranes.

Trajectory-based training enables protein simulations with accurate folding and Boltzmann ensembles in cpu-hours.

An ongoing challenge in protein chemistry is to identify the underlying interaction energies that capture protein dynamics. The traditional trade-off in biomolecular simulation between accuracy and computational efficiency is predicated on the assumption that detailed force fields are typically well-parameterized, obtaining a significant fraction of possible accuracy. We re-examine this trade-off in the more realistic regime in which parameterization is a greater source of error than the level of detail in the force field. To address parameterization of coarse-grained force fields, we use the contrastive divergence technique from machine learning to train from simulations of 450 proteins. In our procedure, the computational efficiency of the model enables high accuracy through the precise tuning of the Boltzmann ensemble. This method is applied to our recently developed Upside model, where the free energy for side chains is rapidly calculated at every time-step, allowing for a smooth energy landscape without steric rattling of the side chains. After this contrastive divergence training, the model is able to de novo fold proteins up to 100 residues on a single core in days. This improved Upside model provides a starting point both for investigation of folding dynamics and as an inexpensive Bayesian prior for protein physics that can be integrated with additional experimental or bioinformatic data.

Accurate calculation of side chain packing and free energy with applications to protein molecular dynamics.

To address the large gap between time scales that can be easily reached by molecular simulations and those required to understand protein dynamics, we present a rapid self-consistent approximation of the side chain free energy at every integration step. In analogy with the adiabatic Born-Oppenheimer approximation for electronic structure, the protein backbone dynamics are simulated as preceding according to the dictates of the free energy of an instantaneously-equilibrated side chain potential. The side chain free energy is computed on the fly, allowing the protein backbone dynamics to traverse a greatly smoothed energetic landscape. This computation results in extremely rapid equilibration and sampling of the Boltzmann distribution. Our method, termed Upside, employs a reduced model involving the three backbone atoms, along with the carbonyl oxygen and amide proton, and a single (oriented) side chain bead having multiple locations reflecting the conformational diversity of the side chain's rotameric states. We also introduce a novel, maximum-likelihood method to parameterize the side chain interactions using protein structures. We demonstrate state-of-the-art accuracy for predicting χ1 rotamer states while consuming only milliseconds of CPU time. Our method enables rapidly equilibrating coarse-grained simulations that can nonetheless contain significant molecular detail. We also show that the resulting free energies of the side chains are sufficiently accurate for de novo folding of some proteins.

A simple approach to spectrally resolved fluorescence and bright field microscopy over select regions of interest.

A standard wide field inverted microscope was converted to a spatially selective spectrally resolved microscope through the addition of a polarizing beam splitter, a pair of polarizers, an amplitude-mode liquid crystal-spatial light modulator, and a USB spectrometer. The instrument is capable of simultaneously imaging and acquiring spectra over user defined regions of interest. The microscope can also be operated in a bright-field mode to acquire absorption spectra of micron scale objects. The utility of the instrument is demonstrated on three different samples. First, the instrument is used to resolve three differently labeled fluorescent beads in vitro. Second, the instrument is used to recover time dependent bleaching dynamics that have distinct spectral changes in the cyanobacteria, Synechococcus leopoliensis UTEX 625. Lastly, the technique is used to acquire the absorption spectra of CH3NH3PbBr3 perovskites and measure differences between nanocrystal films and micron scale crystals.

Raman Vibrational Shifts of Small Clusters of Hydrogen Isotopologues.

Raman vibrational shifts of small parahydrogen (pH2), orthodeuterium (oD2), and paratritium (pT2) clusters with respect to the free molecules are calculated by combining a first order perturbation theory approach with Langevin equation Path Integral Ground State (LePIGS) simulations [ J. Phys. Chem. A 2013 , 117 , 7461 ]. Our theoretical predictions are compared to existing cryogenic free jet expansion results for pure (pH2)N clusters [ Phys. Rev. Lett. 2004 , 92 , 223401 ] and to new measurements for (oD2)N clusters reported here. This method has been successfully used before to predict the Raman vibrational shifts of (pH2)N clusters [ J. Chem. Phys. 2014 , 141 , 014310 ]. The 6-D interaction potential of Hinde [ J. Chem. Phys. 2008 , 128 , 154308 ] is reduced to 1-D using the Adiabatic Hindered Rotor approximation to yield effective pair potentials for both molecules being in the ground vibrational state, and for one of them carrying one quantum of vibrational excitation. These reduced 1-D potentials are fitted to a Morse Long Range analytic form for later convenience. Good agreement between experiment and theory is found for the smaller clusters, but significant deviations remain for the larger ones.

First-principles prediction of the Raman shifts in parahydrogen clusters.

We report a first-principles prediction of the Raman shifts of parahydrogen (pH2) clusters of sizes N = 4-19 and 33, based on path integral ground-state simulations with an ab initio potential energy surface. The Raman shifts are calculated, using perturbation theory, as the average of the difference-potential energy surface between the potential energy surfaces for vibrationally excited and ground-state parahydrogen monomers. The radial distribution of the clusters is used as a weight function in this average. Very good overall agreement with experiment [G. Tejeda, J. M. Fernández, S. Montero, D. Blume, and J. P. Toennies, Phys. Rev. Lett. 92, 223401 (2004)] is achieved for p(H2)(2-8,13,33). A number of different pair potentials are employed for the calculation of the radial distribution functions. We find that the Raman shifts are sensitive to slight variations in the radial distribution functions.

HERV-K-specific T cells eliminate diverse HIV-1/2 and SIV primary isolates.

The genetic diversity of HIV-1 represents a major challenge in vaccine development. In this study, we establish a rationale for eliminating HIV-1-infected cells by targeting cellular immune responses against stable human endogenous retroviral (HERV) antigens. HERV DNA sequences in the human genome represent the remnants of ancient infectious retroviruses. We show that the infection of CD4+ T cells with HIV-1 resulted in transcription of the HML-2 lineage of HERV type K [HERV-K(HML-2)] and the expression of Gag and Env proteins. HERV-K(HML-2)-specific CD8+ T cells obtained from HIV-1-infected human subjects responded to HIV-1-infected cells in a Vif-dependent manner in vitro. Consistent with the proposed mode of action, a HERV-K(HML-2)-specific CD8+ T cell clone exhibited comprehensive elimination of cells infected with a panel of globally diverse HIV-1, HIV-2, and SIV isolates in vitro. We identified a second T cell response that exhibited cross-reactivity between homologous HIV-1-Pol and HERV-K(HML-2)-Pol determinants, raising the possibility that homology between HIV-1 and HERVs plays a role in shaping, and perhaps enhancing, the T cell response to HIV-1. This justifies the consideration of HERV-K(HML-2)-specific and cross-reactive T cell responses in the natural control of HIV-1 infection and for exploring HERV-K(HML-2)-targeted HIV-1 vaccines and immunotherapeutics.